Theory of electron Zitterbewegung in graphene probed by femtosecond laser pulses 



o 
o 

(N 



o 



CD 



o 

> 

m 

00 

O 



X 



Tomasz M. Rusin^Q and Wlodek Zawadzkif 
*PTK Centertel Sp. z o.o., ul. Skierniewicka lOA, 01-230 Warsaw, Poland 
^Institute of Physics, Polish Academy of Sciences, Al. Lotnikow 32/46, 02-688 Warsaw, Poland 

(Dated: July 30, 2009) 

We propose an experiment allowing an observation of Zitterbewegung (ZB, trembling motion) of 
electrons in graphene in the presence of a magnetic field. In contrast to the existing theoretical work 
we make no assumptions concerning shape of the electron wave packet. A femtosecond Gaussian 
laser pulse excites electrons from the valence n = —1 Landau level into three other levels, creating 
an oscillating electron wave packet with interband and intraband frequencies. Oscillations of an 
average position of the packet are directly related to the induced dipole moment and oscillations of 
the average packet's acceleration determine emitted electric field. Both quantities can be measured 
experimentally. A broadening of Landau levels is included to make the description of ZB as realistic 
as possible. Criteria of realization of a ZB experiment are discussed. 

PACS numbers: 73.22.-f, 73.63.Fg, 78.67.Ch, 03.65.Pm 



I. INTRODUCTION 

Zitterbewegung (ZB, trembling motion) of relativis- 
tic electrons in a vacuum was predicted nearly 80 years 
ago by Schrodinger Unfortunately, both the spa- 
cial extension of the ZB motion, being of the order of 
Ac = h/mc, and the ZB frequency ojz — 1ra(? jfi are 
far beyond current experimental possibilities. However, 
it was recently shown that, because of an analogy be- 
tween the behavior of relativistic electrons in a vacuum 
and that of electrons in narrow gap semiconductors 0, [3] , 
one can expect the trembling motion of electrons in nar- 
row gap semiconductors having much more advantageous 
characteristics: the frequency — Eg/h and the ampli- 
tude \z = h/rriQU, where Eg is the energy gap, mg is 
the electron effective mass, and u — (i?g/2TOo)^/^ is the 
maximum electron velocity in the two-band energy spec- 
trum It was further shown that the ZB-like motion 
should occur in other two-band situations, both in solids 
0, H, 0, IMtUUiHi^iP'' "l^' 'isl, [13, m and in other 
systems 



M llil, il 111, i24j, J5,, ili . An observation 



of an acoustic analogue of ZB was reported recently in 
sonic crystals J2l\. 

It was predicted some time ago by Lock [28j] that, if 
an electron is represented by a wave packet, the ZB phe- 
nomenon will have a transient character. This was con- 
firmed by very recent calculations which predicted that 
the decay times of ZB are of the order of femtoseconds 
to microseconds depending on the system in question 
^2, 16; 21]. However, it was also shown that the presence 
of an external magnetic field and the resulting Landau 
quantization of the electron spectrum 'stabilizes' the sit- 
uation making the ZB oscillations stationary in time, if 
one neglects the loss of electron energy due to dipole radi- 
ation [isl . [2^ . It is known that an external magnetic field 
does not induce interband electron transitions, so that an 
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interference of electron states corresponding to positive 
and negative energies remains unchanged and an appear- 
ance of interband frequencies remains the signature of ZB 
phenomenon. On the other hand, due to Landau quan- 
tization of the electron and hole energies also intraband 
(cyclotron) excitations appear in the spectrum. 

All the recent theoretical work on ZB assumed that 
initially the electrons are represented b y G aussian wave 
packets @, HI, M, M, M, S IH, IH, [30| . While this 
assumption represents a real progress compared to the 
initial work that had treated electrons as plain waves 
H i, 0, i, i, [H E mi Hi, it is obviously an idealiza- 
tion since it is not quite clear how to prepare an electron 
in this form. It is the purpose of our present work to 
propose and describe an observation of electron ZB in 
semiconductors that can be really carried out. Namely, 
we calculate a reaction of an electron in graphene excited 
by a laser pulse, not assuming anything about initial form 
of the electron wave packet. In our description we take 
into account currently available experimental possibili- 
ties. Also, we include a broadening of Landau levels and 
investigate its effect on the trembling motion. It is our 
hope that this proposition will help to observe for the 
first time this somewhat mysterious effect that is funda- 
mental for both relativistic electrons in a vacuum and 
electrons in narrow gap semiconductors. 

The following conditions should be met for a successful 
observation of ZB: a) The ZB frequency must be in the 
range of currently detectable regimes, i.e. of the order 
of ujz ~ 1 fs^^, and the size of oscillations should be 
of the order of a few A; b) The ZB oscillations should 
be persistent or slowly transient; c) Both positive and 
negative electron energies must be excited with a suffi- 
cient probability; d) To avoid many-electron effects (see 
Refs. j32, ii4j), the wave packet should be created in a 
one-electron regime. A system that in our opinion ful- 
fills the above criteria is p-type monolayer graphene in 
a constant magnetic field. The wave packet should be 
created by an ultra short monocycle or sub-monocycle 
laser pulse. Because of a very wide frequency spectrum 
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FIG. 1: Electron energy levels for monolayer graphene in a 
magnetic field. Proposed position of the Fermi level is indi- 
cated. Arrows show interband (a), intraband (b), and funda- 
mental (c) energies, see text. 



of such a pulse, the resulting wave packet will have both 
positive and negative energies. The electron oscillations 
give rise to a time-dependent dipole moment which will 
be a source of electric field and it will emit or absorb radi- 
ation in the far infrared range. Experimental parameters 
necessary to create the optical wave packet and to detect 
the radiation should be within the current experimental 
possibilities. 

Our paper is organized as follows. In sections II and III 
we calculate the electron reaction to a short laser pulse, 
the creation of ZB, and we describe an electric field of 
radiation emitted by the trembling electron. In Section 
IV the influence of Landau level broadening on the trem- 
bling motion is investigated. In Section V we describe 
the time-dependent luminescence filtered by a time gate 
and a frequency filter. Finally, we discuss our results and 
conclude by a summary. 



II. PRELIMINARIES 

We consider p-type monolayer graphene in the pres- 
ence of a magnetic field perpendicular to the layer. The 
electron- hole Hamiltonian Hi at the Ki point of the Bril- 
louin zone is given by [s^ [s^ [sj 



Hi = -Ml 



a 



(1) 



where 



are the lowering and rising operators and 



the characteristic frequency of the system is 57 = y/2u/ 
where L = y^h/eB is the magnetic length, and the ve- 



locity u ~ 10^ cm/s. The energy spectrum of Hi is 
En ~ sgn(n)/iO-\/|ri[, where n — 0, ±1,..., see Fig. 
[1] The eigenstates of Hamiltonian ([1]) for the gauge 
A = [-By, 0,0] are 



-sgn(n)(/)|„|_i(C) 
4>\n\{i) 



(2) 



where ^ = y/L-fc^jL and 0„(^)=(C„/V^) e 2C^H„(^) is 
the n-th eigenstate of the harmonic oscillator. Here C„ = 
l/y^^hTJy^ , and H„(f) are the Hermite polynomials. 
For n = 0, the first component in Eq. ^ vanishes and 
the normalization coefficient is 1/ \p2j:. We assume the 
Fermi level to coincide with the Landau level (LL) n = 
—2 and consider the initial electron in n = — 1 state, see 
the Discussion below. 

The wavelength of the laser light is assumed to be much 
larger than the spacial size L of the n = — 1 state, so we 
can neglect spacial variation of the electric field in the 
laser pulse. We take the perturbing potential due laser 
light in the form 



W{t) = -ey£oe-^' cos(a;Li), 



(3) 



where e is the electron charge, r is the pulse duration 
(FWHM), ul — 27rc/AL is the laser frequency (being of 
the order of 3xl0^^s~^), and £o is the amplitude of elec- 
tric field. A Gaussian shape of the laser pulse is widely 
used in optical experiments and it parameterizes effec- 
tively a profile of electric field in the laser beam. 

As a result of a laser shot, the initial state of the 
system ^k{t) — i/'fee"'^''*^^ evolves into the final state 
^fe(0 — X]j '^j(OV'fce~*"^^*^'', which is a combination of 
the eigenstates of Hi with suitably chosen coefficients 
Cj{t). The resulting time-dependent dipole moment is 
D{t) = e{-^u{t)\r\^k{t)). 

The total Hamiltonian, including the perturbation due 
to the laser light, is 



H = Hi + W{t). 



(4) 



The corresponding time-dependent wave functions are 
^k{t) = e"'-'^*/''*fe(0), and the dipole moment is 



= e(*,(0)|r(<)|vE',(0)) =e(r(t)). (5) 



D{t) = e(*fc(0)e*"*/'*|r|e-'^*/"*fc(0)) 



Here r{t) is the electron position in the Heisenberg pic- 
ture. Thus the dipole moment D(t) is proportional to 
the time-dependent position averaged over the electron 
wave packet. 

A time-dependent dipole moment is a source of elec- 
tromagnetic radiation. We treat the radiation classically 
38 1 and take the radiated transverse electric field to be 
38| 



D{t) s\n{e) 
47reoc2 R ' 



(6) 



3 



where eo is the vacuum permittivity, is an angle be- 
tween the direction of electron motion and a position of 
the observer R. Since D[t) = e{r{t)), Eq. ([6]) relates the 
electric field of the dipole with the average acceleration 
of the packet. If the electric field is measured directly by 
an antenna, one measures the trembling motion of the 
wave packet. If the square of electric field is measured in 
emission or absorption experiments, the signature of ZB 
is the existence of peaks corresponding to interband and 
intraband frequencies and their dependence on packet's 
parameters. Accordingly, in the time resolved lumines- 
cence experiments it should be possible to detect directly 
the motion of the packet with interband and intraband 
frequencies. 



III. EMITTED ELECTRIC FIELD 

Now we calculate explicitly the electric field emitted by 
the trembling electron. If we consider the perturbation 
of Eq. ([4]), the standard time-de pen dent perturbation 
theory gives for the wave function [40] 

^^(t) i,^e-^^-'l^ + ^ cl^(t)i,,e-^^^'l^ + . . . , (7) 

3 

with 



Setting the initial state to be n = — 1 and using Eq. ([3]) 
we have 



1 / V^eL^oT \ y 



where 



a^y = / vLi(r)#,(r)dr, (10) 



-at'Vr%*a;-i,t' cos{uj Lt')dt' . (11) 



Here cu-ij = {E^i - Ej)/h and a = 21n2 « 1.34. Since 
the size of n ~ —1 state is of the order of L, the coef- 
ficients a^_^j are of the order of unity. It is also is easy 
to show that < 1. The dimensionless perturbation 

expansion parameter is C = ^/TreLSoT/hy/a. 

To obtain the coefficients a_ij from Eq. pO|) we cal- 
culate the matrix elements of y between different eigen- 
states ipnir) of Hi. The selection rules for {ipn\y\'ipj) 
of Eq. PH)) are \n\ — \j\ = ±1, so for rt = — 1 there 
are three non-vanishing matrix elements corresponding 
to j = 0, ±2, see Fig. [TJ The approximate wave function 
*_!(<) is then 



-iE-it/h 



i=0,±2 



g-»_E,t/B^^.. (12) 



Therefore, the laser shot creates a non-stationary wave 
packet given by Eq. (HH). This wave packet contains 
states with positive and negative energies, which is a 
necessary condition to create the ZB motion [l^] ■ If 
the packet contains only positive or only negative energy 
states, the ZB will not occur. To calculate D{t) we aver- 
age y and a; over the wave function (jl2p . As a result we 
obtain 16 terms, of which one term does not depend on 
c_ij , six terms are proportional to c_ij, and the remain- 
ing nine terms are of the second order in C-ij. Since the 
zero order term does not depend on time, we concentrate 
on the time-dependent terms of the lowest order in Sq, 
and we have 



1, all, e 



h.c. 



Dy{t) ~ const + eL c\ 

j=0,±2 

(13) 

and similarly for Dx(t), with a^ij replaced by a^ij- Be- 
cause the pulse duration r is much shorter than the 
period of ZB oscillations Tz ~ 2tt/Q, we approximate 
b-ij{t) by b-ij{oo). Then [see Eq. (HH)] 



^ s=±\ 

Within this approximation we have 



(14) 



-10 . 



sin(a;gt) + B 6_i2 sin(u;f i) 



-B+ fe_i„2 sin(tjjt)) , 



5-1 



= rfo ( -y^ cos(wot) + B fe_i2 cos(tjf i) 
S+ 6_i_2Cos(cjJt)) , 



(15) 



where = fi(Vrr+T — ^/n), — Q.{\Jn + 1 + -/tT), 
B^ = \/t/2 ± 3/4, and do = -eLC,. Taking the second 
time derivative of the dipole moment we find the electric 
field components of the emitted electromagnetic wave 



f,(r, t) = S(r) ( ^ sin(c.gt) -f ^ sin(c.f i) 



2 



sin(a;Jt) 



4 

E^{r, t) = S(r) ( - cos(wgt) + cos{uj^t) 



b-i- 



■ COs{u!it) 



(16) 



where S(r) = do^l^ sin(6')/(47reoc2i?^). Equations (fTS]) 
and (I16p are among the main results of our work. They 
state that both the induced dipole moment and the cor- 
responding electric field oscillate with three frequencies. 
The frequency w_i2 = {V2+ 1) ft corresponds to the Zit- 
terbewegung, i.e., to the motion of the packet with an in- 
terband frequency. This frequency corresponds to the in- 
terband ZB frequency ujz = 2mf,c^ /h of relativistic elec- 
trons in a vacuum. The interband frequency is character- 
istic of ZB because the trembling motion occurs due to 
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FIG. 2: Oscillations of dipole moment during the first f 000 fs 
of electron motion after the laser pulse. Experimental charac- 



W/cm^ 



r = 1.6 fs, B = 1 



teristics: pulse intensity Ix 10 
T, b) r = 1.6 fs, B = 10 T, c) r = 3.0 fs, B = 1 T, d) r = 3.0 
fs, B = 10 T. Dipole moments in a) and b) are in 10"^** [Cm] 
units, while in c) and d) they are in IQ"^'^ [Cm] units. The 
above results refer to very narrow Landau levels, disregarding 
broadening due to electron scattering and the e-e interaction, 
see Section IV. 



an interference of electron states related to positive and 
negative electron energies [3l|, [s^l • The second frequency 
<j;_i_2 = (\/2 — 1) ^ describes the intraband cyclotron 
motion of the packet. The third frequency to-io = 
has both interband and intraband character (see Fig. [1]). 
Contrary to the relativistic quantum mechanics, in zero- 
gap materials like graphene the interband ZB frequency 
is not vastly larger than the cyclotron frequency. 

In Fig. [5] we plot the oscillating dipole moment within 
the first 1000 fs of motion after the laser shot for two 
magnetic fields B — 1 T and B = 10 T, and two laser 
pulses. The first pulse (Figs. [5^, O)) has a duration 
T — 1.6 fs and a base laser wavelength = 650 nm. 
This pulse has a sub-monocycle duration and it is the 
shortest pulse created experimentally within the visible 
laser wavelength 4]J. In Figs. O;, we assume pulse 
duration t — 3.0 fs and a laser wavelength Xl — 720 nm. 
This pulse has 1.25 of the laser monocycle and its exper- 
imental properties were discussed in Ref. [12]. The use 
of a few monocycle pulses (with r > 5 fs) is not effective, 
since the probability of excitation of a wave packet in Eq. 
(jl2p drops exponentially with pulse duration r, see Eq. 

In Fig. [3] we plot the corresponding electric field for 
the same parameters during the first 250 fs of oscilla- 
tions. We assume the laser intensity to be / = l.Ox 
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FIG. 3: Electric field emitted by one electron during the first 
250 fs of electron motion after the laser pulse. Experimental 
characteristics: pulse intensity Ix 10^ W/cm^, a) r = 1.6 fs, 
B = 1 T, b) r = 1.6 fs, B = 10 T, c) r = 3.0 fs, B = 1 T, 
d) r = 3.0 fs, B = 10 T. Note the difi'erence in magnitudes 
for T = 1.6 fs and r = 3.0 fs. The above results refer to 
very narrow Landau levels, disregarding broadening due to 
electron scattering and the e-e interaction, see Section IV. 



10^ W/cm^, the emitted electromagnetic wave detected 
at the angle 9 = 45°, and the distance R — 1 cm. All 
the quantities in Figs. [2] and [3] are calculated per one 
electron. Since the frequencies are incommensurable, the 
electron trajectories r(i) are not closed and there is no 
repeated pattern of oscillations. The motion of the wave 
packet is permanent in the time scale of femtoseconds or 
picoseconds but there is damping of the motion due to 
the light emission in a long time scale. The results shown 
in Figs. [2] and [3] refer to very narrow Landau levels, dis- 
regarding broadening due to electron scattering and the 
e-e interaction, see Section IV. 

We can draw the following qualitative conclusions from 
Figs. [2 and [21 First, for small magnetic fields B the pe- 
riod of oscillations is longer than for large fields, which 
is related to the basic frequency fl. Second, irrespective 
of the variation of with B, for small fields the oscil- 
lations are dominated by the low (cyclotron) frequency, 
while at stronger B the high (ZB) frequency dominates. 
Finally, comparing the magnitudes of dipole moment or 
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TABLE I: Coefficients h-xj for the electric field in Eq. (HH) 
for different pulse durations r and various magnetic fields B. 



which corresponds to /a = VI/{2tt) = 8.61 THz. The 
cyclotron frequency is fc — 3.53 THz, while the ZB fre- 
quency is fz = 20.8 THz. All the three frequencies are 
within the range of currently available THz photocon- 
ductive antennas, see e.g. Ref. [4^], and it should be 
possible to detect the emitted field experimentally. In 
contrast, for _B = 40 T the corresponding frequencies are 
fa = 54.5 THz, = 22.3 THz and fz = 131 THz, and 
they are more difficult to detect. This is the other reason 
for using low magnetic fields in the experiment. 



emitted electric field for t = 1.6 fs with the correspond- 
ing values for r = 3.0 fs we observe that the amplitude 
of oscillations depends very strongly on the duration r of 
the pulse. 

To analyze these effects quantitatively we collected in 
Table 1 the coefficients h-ij of Eq. (fT4|) . used for the 
calculation of electric field in Eq. (|16p . The results pre- 
sented in Table 1 show that for fixed B and r = 1.6 fs 
the coefficients 6_i_2 and 6_io are nearly three orders of 
magnitude larger than those for r = 3.0 fs. We also note 
that, for i? = 1 T, all are nearly identical, while 
for B = 40 T there are visible differences between vari- 
ous h-ij. This difference is a factor of 3 for t = 1.6 fs, 
while for r ~ 3.0 fs the coefficient 6-21 is two orders of 
magnitude larger than 6-1-2- This explains the domi- 
nance of the interband ZB frequency for large B in Figs. 
[5] and [3] The conclusion from this analysis is that the 
optimum conditions for observing the ZB, i.e. the packet 
motion with both interband and intraband frequencies, is 
the regime of magnetic fields of a few Tesla, since in this 
regime the two kinds of motion exist with comparable 
weights. One should note that the coefficients as 
defined in Eq. (fTTjl , are closely related to the coefficients 
c-ij in the perturbation series, see Eqs. ([9]) and (fT4|) . A 
practical lower limit for magnetic fields is the condition 
that an energy distance between LLs should be larger 
than the widths Tl of LLs in graphene. According to 
Ref. [43] one observes resonant magneto-optical transi- 
tions, both interband and intraband, in graphene begin- 
ning with a magnetic field of B k, 0.4 T. Thus B « 0.5 
T seems to be the lowest possible magnetic field suitable 
for the experiment described above. The Landau level 
broadening and its effect on ZB are discussed in the next 
section. 

Comparing Figs. [2] and [3] for fixed r one observes a 
strong dependence of the dipole moment on the magnetic 
field, while the radiated electric field depends only weakly 
on B. The reason is that the electric field of Eq. (fTB]) is 
multiplied by a factor cIqVL^ . Because do ~sLC., C, <x L, 
and n = y/2u/L, the product daVL^ does not depend on 
B. On the other hand, do oc oc B~^, and the dipole 
moment D{t) depends strongly on B. Therefore, the 
electric field of the emitted electromagnetic wave (and 
the radiated power) does not change significantly with 
the magnetic field intensity, it depends on B only via 
coefficients see Table 1. 

For i? = 1 T, the basic frequency fl is 5.41x 10^^ s~^. 



IV. ZB IN REAL SAMPLES 

In the previous section we considered an idealized case 
of very narrow Landau levels in graphene. In real sam- 
ples additional effects occur and their presence affects the 
motion of the wave packet. Two effects may play a role 
in the proposed experiment: the electron - electron (e- 
e) interaction ^ and the presence of disorder [i^, [ivl . 
The scanning tunnelling spectroscopy results of Ref. [43 
indicate that the e-e interaction leads to a Lorentzian 
shape of DOS of the Landau levels and it opens an energy 
gap between the electrons and holes. Thus the massless 
Dirac fermions acquire a small non-zero mass. As shown 
in the numerical simulations of Ref. [4^, the presence 
of disorder changes the shape of DOS from Lorentzian to 
Gaussian peaks. Additionally, the disorder potential may 
change the position of the Fermi level within sample. 

The band-gap caused by the e-e interaction is of the 
order of 10 meV [1^, so at a magnetic field of i? = 1 
T it is much smaller than the basic energy hV, « 36 
meV. In this case the energy spectrum of graphene in 
a magnetic field is described by an analogue of the Dirac 
equation, whose energy levels and eigen-functions are 
well known. The trembling motion of the packet will 
not change qualitatively, as compared to the above de- 
scription, but it will oscillate with the interband fre- 
quency Cl = y^fi^ _|_ Eg^p/h?, and all frequencies will 

have slightly different values than those calculated in the 
gapless model. 

On the other hand, the broadening of the Landau levels 
may strongly influence the trembling motion of the wave 
packet. To analyze the overall impact of all the effects 
leading to the level broadening: disorder, e-e interaction, 
electron-phonon scattering, etc., we assume finite widths 
of all energy levels, characterized by broadening param- 
eters r„. We treat r„ as phenomenological quantities 
determined experimentally and including all scattering 
mechanisms existing in real samples. We approximate 
the broadening of DOS by a Lorentzian line-shape (isj 
irrespective of the detailed scattering mechanism [46| . In 
this approximation the Landau energies En are replaced 
by complex energies En = En + ir„. After the replace- 
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FIG. 4: Calculated electric fields Ex{t) and Ey{t) emitted by 
one trembling electron during the first 400 fs after the laser 
pulse. Pulse parameters: intensity Ix 10^ W/cm^, r = 1.6 
fs, magnetic field B = 1 T. Bold lines - electric fields for 
broadened Landau levels described by Lorentzian line-shapes 
with experimental values of r„ (see text). Thinner lines - 
electric fields for delta-like Landau levels. The latter results 
are the same as those shown in Fig. 



ment the dipole moment of Eq. (fT3|) changes to 
Dlit) ~ const + eL ^ cL„ a^y e''"-!^*"^^* +h.c. + . . . , 

(17) 

which leads to [cf. Eq. ^5)) ] 

Dlit) ^do ^sin(^St)e-i^«* 

z.\-r.nt 



+B b-i2Sm{ujft)e 
+B+ 6-1-2 sin(wit)e 



Dlit) =do (^cos(^„^t)e-i^°* 

B- fe-i2Cos(wft)e-^+^* 
B+ 6-1-2 cos(wJt)e-^-2*) 



(18) 



As to the numerical values of r„, we take after Ref. [45[ 
r-i-2 = 5.1 meV. For the state n = we use Tq = 4.3 meV 
pO | which corresponds to the disorder potential Vg = 120 
meV [13]. These values are similar to the line widths of 
r„ = 7 meV measured at i? = 1 T in far infrared trans- 
mission experiments 48|, and r„ = 1.6 meV determined 
by the Quantum Hall Effect at B - 40 T [11]. 

The electric field E{t) is calculated, as before, as a sec- 
ond time derivative of D^{t), see Eq. (|15p . In Fig. [3]we 
plot the electric field emitted by an oscillating electron 



FIG. 5: Intensity of three emission lines versus magnetic field 
for two laser pulses of different durations t. 



within the first 400 fs of motion after the laser shot of 
the width t — 1.6 fs in a magnetic field B = 1 T. The 
two bold lines describe calculated electric fields Ex{t) and 
Ey (t) for the Landau levels having the broadening param- 
eters r„ indicated above. The two thinner lines show the 
electric fields calculated without damping (r„ = 0), see 
Eq. dTS]) and Fig. [3ji. Within the first 50 fs of motion 
the electric fields emitted in the two cases are similar, but 
later the damping of the emitted fields for broadened lev- 
els is visible. After around 400 fs the trembling motion in 
real case disappears. It can be seen that the maxima of 
oscillations for the damped ZB motion coincide with the 
undamped ones. The general conclusion from Fig. 2] is 
that the existence of disorder, many-body effects or other 
scattering mechanisms changes the persistent ZB motion 
to a decaying one, within the characteristic lifetimes for 
these processes: r„ = l/r„ « 130 fs. Nevertheless, since 
the parameters F^ used in the calculations correspond to 
the measured lifetimes in real graphene samples, it fol- 
lows that the broadening of the Landau levels does not 
prevent the existence of ZB. Clearly, a lower disorder in 
better samples will result in longer decay times for ZB. 



V. TIME-RESOLVED LUMINESCENCE 

Knowing the electric field emitted by an oscillating 
dipole we can calculate the intensity of radiated light. 
If £{t) is given by a sum of cosine (or sine) functions: 
£{t) = fj cos{LUjt), the emitted power averaged over 

a long time is P = (1/2) J2j l/jP- The total power S 
passing through a closed spherical surface at a distance 
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CO (ps" ) 



FIG. 6: Calculated power spectra of the emitted light for 
different opening gate times T. For the assumed resolutions 
7 = 3 cm~^ and F = 1000 cm~^ the peaks At uj — Q, and 
Lo — (\/2~l) ^ are not resolved. Inset: Intensity of two peaks; 
Lu = Q. and — (\/2 + 1)0, versus gate time T. Repetition 
of oscillation pattern occurs after Ta ^ V2 Q. ^ 82.2 fs for 
B = 1 T. 



R from the sample can be calculated integrating P over 
the enclosing surface. Using Eq. (|16p and the above 
formula we find 



S = 



(4|6_ 



101 



121 



-1-21 



(19) 



which is the Larmor formula for our problem. The emit- 
ted spectrum S consists of three lines of different inten- 
sities having the frequencies lo — uj — (\/2 + 1)0 
and LO = (-^72 — 1)0,. The existence of lines with inter- 
band frequencies is the signature of Zitterbewegung in 
this system. In Fig. [S] we plot the relative light inten- 
sities for different emission lines versus magnetic field B 
on the logarithmic scale. The upper curves correspond 
to r = 1.6 fs, the lower ones to r = 3.0 fs. In both cases 
the intensities depend on magnetic field, this dependence 
is most pronounced for the ZB frequencies. However, for 
r = 1.6 fs the dependence of intensities on the magnetic 
field is much weaker than for r — 3.0 fs. At high B the 
spectrum is dominated by the ZB frequencies. It should 
be noted that the intensity of radiation is proportional 
to O** and to cx £q, as for the Thompson scattering. 

To observe the motion of electron represented by a 
wave packet one can use the time-resolved luminescence 
[sol [Sll . [5^ . In this technique, the electric field E{i) of 
the light emitted by a sample is transmitted through two 
filters: a time gate B{t,T) 



and a frequency filter H{t,ujp) 



H{lo,lof) = 



"f^ + {lo ~ ufY 



(21) 



The time gate lets the field pass for T — r/2 < t < 
T + r/2, where F is the gate's parameter describing the 
width of the window, and T is the center of the window. 
The frequency filter selects frequencies close to LOp with 
the resolution 7. The transmitted field is 

SHB{t) = I H{t-t',LOF)B{t',T)£{t')dt', (22) 



where H{t^ujp) is the Fourier transform of H{lo,lof). 
The time-dependent spectrum of £{t) is 



S{T,wf) = fdA f \£HB{t)\^ dt, 



(23) 



where the first integration is over the sphere of radius 
R enclosing the sample. For £{t) given by Eq. we 
obtain 



S{T,lof) = M 



'^^)^G(«,/3), (24) 



where M is a constant and a,l3 — {n,j,s). The sum- 
mation is restricted to n = —1, j — 0,±2, and s = ±1. 
The index r = x,y, and we denote uja,ijji3 = swnj- Fac- 
tors and bp are the coefficients in front of trigono- 
metric functions in Eq. (fTB|) . We have h^^r 



12 



6_i2/4, and 



-&-1-2/4. 



B(t,r)=exp(-F|t-r|), 



(20) 



- 6-10/2, 
Similarly 

&iio = -b-io/2, 6^12 = ^'-12/4, and b^_^_,^ = 6_i_2/4. 
The analytical function G{a, 0) is defined by Eq. (17) of 
Ref. 

In Fig. [S]we show the calculated time-dependent spec- 
trum of ZB oscillations for the following parameters: 
laser pulse width t = 1.6 fs, laser wavelength Xl = 650 
nm, magnetic field _B = 1 T, frequency filter 7 = 3 cm~^, 
and the gate width F = 1000 cm^^. Three curves cor- 
respond to frequency spectrum observed at three gate 
opening times T. If follows from Eq. p4)) and from the 
selection rules for j, see Eq. (|10p . that the time depen- 
dent spectrum S{T,lof) has maxima at the frequencies: 
LO — LOc, LO — f2, and lo = loz- Because of finite res- 
olutions of the time gate and the frequency filter, the 
maxima for lo = lOc and lo = form one unresolved 
peak. 

Changing the gate opening time T one should ob- 
serve different frequency spectra. The spectrum re- 
peats its pattern after Ta ~ 27r/cja = 82.2 fs, where 
LOa = LO-io — oJ-1-2 = V^Q. To iUustratc this time de- 
pendence we plot in the inset of Fig. [5] intensities of 
the emitted radiation for two frequencies: lo — and 
LO — VLz versus gate's opening time T. We see that, for 
both maxima, the intensities oscillate with the period T^. 
The oscillating pattern of S{T,u!f) is also the signature 
of Zitterbewegung. 
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VI. DISCUSSION 

In section II we assumed that the Fermi level in mono- 
layer graphene coincides with the Landau level n = — 2 
and considered the initial electron in the state n = — 1, 
see Fig. [TJ This means that, before the electron re- 
acts to a short laser pulse, it must be pumped to the 
state n = —1 from lower Landau levels. A conventional 
light source is suitable for such pumping but, in order to 
achieve high intensities of the emitted lines, one should 
use a laser pump in resonance with huj = E_i — E_2 en- 
ergy. It should be emphasized that the upper component 
of the state n = — 1 in a magnetic field is described by 
the Gaussian wave function in space, but, since it is an 
eigenstate, it does not have a time dependence. The de- 
cisive factor is the subsequent laser pulse which excites a 
series of electron eigenstates. 

The results described by Eq. (fT6|) are obtained for one 
electron. The population of the Landau level is eB/h 
and the total intensity of radiation is obtained by adding 
the contributions from all electrons in the initial n = — 1 
level. In our proposition we select the initial electron 
state n = —1 for several reasons. Due to the selection 
rules A|j| = ±1 there exist three frequencies Lo^ij con- 
tributing to the electron motion. If the state n — Q were 
selected, there would be only two non-vanishing matrix 
elements yo,-i = yo,+i corresponding to the same fre- 
quency ri, see Fig. [TJ Assuming the Fermi energy at 
the n = —2 LL and supposing, in consequence, the state 
n = —2 to be roughly half filled, we avoid the Pauli ex- 
clusion principle in the calculation of the dipole matrix 
elements. Therefore, one-particle formalism of ZB cal- 
culation can be applied, see [H, [S^l- The presence of 
a magnetic field is essential for a successful experiment 
since, as we said above, in the absence of magnetic field 
the ZB oscillations have a transient character with the 
decay time of tenths of femtoseconds [l^, [H, [3 and the 
detection of such oscillations is difficult. However, as we 
showed in Section IV, the broadening of Landau levels 
also leads to a transient character of ZB. In high qual- 
ity graphene samples with small disorder the electron ZB 
should last longer after the laser pulse. 

From the early eighties, short laser pulses were used in 
quantum chemistry to excite wave packets in molecules 
[53]. After a laser shot, the non-stationary wave packet 
evolves in time and its motion is measured in many ways: 
absorption, luminescence, Raman scattering, etc., using 
the so called pump and probe method. The experiment 
proposed here is in principle similar to these techniques. 
Our additional requirement is the necessity to use mono- 
cycle or sub-monocycle laser pulses in order to excite the 
electron packet with both positive and negative energies 

In the nineties, the Bloch oscillations were observed 
in superlattices in a static electric field. In these exper- 
iments the electric field creates a set of discrete levels 
(Wannier-Stark ladder, WSL). Then, a short laser pulse 
creates a wave packet consisting of many discrete states 



and an oscillating dipole moment appears. The oscil- 
lations of the dipole moment were measured either by 
detection of THz radiation |54| , or by measuring change 



in the positions of WSL energies in pump and probe ex- 
periments [55]. There exist several similarities between 
the Bloch oscillator experiment and our proposition. In 
both cases the system is quantized by an external field, 
the laser pulse creates an electron wave packet consisting 
of many discrete levels, and the oscillations of the wave 
packet lead to the time-dependent dipole moment that 
can be observed experimentally. The main difference is, 
again, that in order to observe ZB one needs to create an 
electron packet having both positive and negative energy 
states. Finally, similar techniques to the one proposed 
here were applied to observe a coherent emission from 
double- well systems [H, 113] ■ 

In our treatment we used the first order time- 
dependent perturbation theory, see Eq. (|12p . For very 
strong electric fields this approach may be insufficient. 
However, the fast decrease of coefficients with in- 
creasing pulse width r assures in practice the validity of 
the perturbation expansion. 

It was shown in the previous work on ZB that its exis- 
tence is related to a nonzero momentum of the electron 
[l,[l2]. In our gauge for A, the initial electron state has 
an initial momentum k^- A laser pulse W{t) creates the 
state with nonzero momentum also in the y direction. 
Direct calculations show that at t = the average mo- 
mentum is {pi)= {'^{Q)\{h/i){d/dxi)\'^{Q)) which gives 
{Py) ~ '2'C,fi/ L and (px) — k^. These results are obtained 
in the lowest order in £q. The asymmetry between x 
and y directions follows from the electric field E directed 
parallel to the y axis. 

As we said above, several conditions should be met in 
order to observe the ZB. In our understanding, all these 
conditions can be fulfilled: p-dopcd monolayer graphene, 
ultrashort pulse of the required intensity, detection of 
the emitted electric field in the 3-21 THz range, or 
an observation of the time resolved spectra. It should be 
possible to satisfy these conditions in a single experiment. 



VII. SUMMARY 

To summarize, we proposed and described a possible 
method to observe the trembling motion of electrons in 
graphene in a magnetic field. The central point is that 
we did not assume anything about the shape of the initial 
electron wave packet. We calculated the time dependent 
dipole moment induced by an ultra-short laser pulse. For 
electrons located initially in the n — ~1 state the in- 
duced dipole moment oscillates with three frequencies, 
of which the frequency uiz = (v^ -I- 1) is the signature 
of Zitterbewegung. A possibility of performing such an 
experiment and detecting the ujz frequency are discussed 
and it appears that the current experimental techniques 
are sufficient for a successful observation of ZB in high 
quality graphene samples. 
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APPENDIX A 

It is well known that in the Brillouin zone of monolayer 
graphene there exist two inequivalent minimum points 
Ki and K2 [13, HI] ■ Above we calculated contributions 
from an electron at the point Ki defined by Hamiltonian 
(H]). The question arises: what is the contribution of 
an electron at the point 7^2? To elucidate this issue we 
consider the Hamiltonian for electrons at the K2 point of 
BZ ^,[11 

H,^nn^l';y (Ai) 

Comparing Eq. (jAip with Eq. ([4]) it is seen that H2 — 
—Hi ■ The eigenenergies of both Hamiltonians are the 
same, but the eigenstates of H2 are 

6**^="^ / (f)\n\{0 \ 

^"^""^^ Vl^lsgnWJ>H-i(e) j' ^^^^ 

which differs from tp{r) of Eq. ([2]) by an exchange of 
the upper and lower components and by the change of 
sign of state. The matrix elements are {iplylip) 

— {x\y\x)^ so both Hamiltonians give the same selection 



rules. The perturbed wave function T(r,t) for the elec- 
tron at the K2 point is [see Eq. p^ ] 

T_i(0 e-'^^-*/Vi + cii, e-^^^-^/'^x,, (A3) 

j=0,±2 

with the same coefficients C-ij as in Eq. (fT2l) . In conse- 
quence, the wave packet at t = for the electron at the 
Ki point is [see Eq. p^ ] 

vi/_i(0) ~ + 

while for the K2 point the initial wave packet is 

Thus, the same laser shot creates two different electron 
wave packets at Ki and K2 points. This is in contrast 
to the assumption made in Ref. [l5| . where the same 
wave packet was assumed for both Ki and K2 points. 
The result of this assumption was a partial cancella- 
tion of one of the electric current components. This 
was unphysical since it violated the rotational symme- 
try of the X — y graphene plane. In the present ap- 
proach, the two wave packets evolve according to different 
Hamihonians. There is -i{t) = e"'^i*/'"'1'_i(0) and 

T_i(f) = e-*^2*/'^T_i(0). A direct calculation shows 
that now the contributions to the electric current and 
dipole oscillations arising from electrons excited at the 
two nonequivalent points Ki and K2 are equal. 
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